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Abstract 

£> | We derive a reation between four-fermion QED Green functions of different co- 

variant gauges which defines the gauge dependence completely. We use the derived 
£T) • gauge dependence to check the gauge invariance of atom-like bound state calculations. 

CN | We find that the existing QED procedure does not provide gauge invariant binding 

' energies. A way to a corrected gauge invariant procedure is pointed out. 



C^; 1 Introduction 

The persisting discrepancy between theory and experiment for positronium width [[jj] is a 
chalenge for QED. At the moment the hope is on taking into account corrections of relative 
order a 2 P, ||. In the circumstances the question of self-consistency of the calculations, in 
particular, of gauge invariance of the result is of prime concern. 

The modern way to calculate parameters of two-particle atom-like bound states is to 
extract them from corresponding four-fermion QED Green function (see, for example, @, |5|, |(| 
and this paper below). Thus, to check the gauge independence of the calculated bound state 
parameters, one should carry the gauge parameter through all the extraction procedure. (An 
example of this see in |7| where the gauge independence of the correction to the positronium 
width of relative order a was checked.) The extraction procedure gets more and more 
complicated with an increase in order of radiative corrections and direct order by order check 
of gauge invariance becomes impractical as a check of self-consistency of the calculations. 
Instead, one would like to exploit gauge invariance choosing a most convenient gauge and 
switching from one gauge to another in the process of the calculations. In view of these 
complications, it seems pertinent to make a step out of the concrete practice of bound- 
state calculations and to study first the gauge dependence of the four-fermion QED Green 
function itself without taking into account the complications of the bound-state parameter 
calculations. 
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In the present paper we derive a relation between four-fermion QED Green functions 
of different values of gauge-fixing parameter (we consider the covariant gauges only). The 
relation completely defines the evolution of the Green function in the gauge-fixing parameter. 
Our derivation does not use perturbation theory. Next, we use our relation to check gauge 
invariance of the extraction procedure of atom-like bound-state parameters. The result 
is negative. It turns out that the existing procedure provides gauge-dependent answers 
for binding energies. We find a flaw in the procedure which is responsible for the gauge- 
dependence of the result and point the way to its correction. 

Next section contains a derivation of the evolution in the gauge-fixing parameter; section 3 
comprises a brief recall of the extraction procedure and an utilization of the general evolution 
formula from section 2 for an analysis of gauge- dependence of the extraction; in the last, 
fourth, section we point out the reason for the gauge dependence and the way to the correct 
procedure. 



2 Evolution in Gauge-Fixing Parameter 

Let us consider the four-fermion QED Green function 

Gp(xf,Tf,Xi,Xi) = i I DipDA exp(iS QE D((3)) (■ip(xf)ip(x f ))(ilj(xi)ilj(xi)) , (1) 



where Xf(Wf) is a coordinate of outgoing particle (antiparticle) and Xj(xj) is the same for 
ingoing pair. The definition of gauge fixing parameter (3 is given by corresponding photon 
propagator: 



dk ( ^k^k v \ i ikx 



/dk I 
7^4 l-fa+P 



(2tt) 4 \ k 2 J k 2 

Our aim is to study the dependence of Gp on (3. To this end, it is useful to consider 
a Green function in external photon field, G(A), which is a result of integration over the 
fermion field in the rhs of (|Tj). From the one hand, it is simply connected to the Green 
function ||: 

Gp = (e L ?G(A)) A=0 , Lp = J£- ( 3 ) 

(In this formula each Lp generates a photon propagator; the dependence on the coordinates 
of ingoing and outgoing particles is suppressed for brevity.) From the other hand, G(A) is 
siply connected to a gauge invariant object Gi nv (A): 

G(A) = G inv {A) exp lie [_ ' A^dx^ - ie f * A^dx^ ) . (4) 

The gauge invariance of Gi nv means that it is independent of the longitudinal component of 
A: 

d^—G mv (A) = (5) 
and is a consequence of gauge invariance of the combination 

ij(x)exp(iej A^Jipiy). (6) 
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A substitution of (Q) into ([3]) yields 

G/3 = I e L "G inv {A) exp lie H A^dx* - ie H A^dxA ) . (7) 

V V Jx f Jx * / ) A=0 

Let us take a (3- derivative of both sides of this equation: 



dp 



G/3 = (e L ^{d p L fi )G inv {A) exp [ie [_ 1 A^dx* -ie T* A^dx" 

\ \ JXf J Xi , 



A=0 

To get an evolution equation, one needs to express the rhs of this equation in terms of Gp. 
It is possible because (dpLp) commutes with G inv (A) and gives a c-factor when acts on the 
consequent exponential. So, (§) transforms itself into 

where we have restored the x-dependence and used F to denote the action of (dpLp) on the 
exponential: 

(dpLp)exp (ie f * A^dx^ — ie f A^dx^ = 

F(xf, Xf, Xi, Ti) exp ( ie j A^dx^ — ie j A^dx^j . (10) 

An explanation is in order: In deriving (§) we have used a commutativity of (dpLp) and 
Gi nv (A); it is a direct consequence of gauge invariance of Gi nv (see (|J)) and the fact that 
(dpLp) contains only derivatives in longitudinal components of A (see (^|) for a definition of 
L/3 and (0) for /5-dependence of D^ v ). 

The solution of eq.(|9|) for /3-evolution is 

Gp(xf,Xf,Xi,Xi) = exj)((^ - p )F(xf,Xf,Xi,x i ))G l 3 (xf,Xf,Xi,Xi). (11) 

To get the final answer one needs an explicite view of F from (pTT|) . It is easily deduced 
from the F-definition (10) and the following representation for the longitudinal part of the 
photon propagator: 

d p D^(/3,x) = --^d ll d l/ \n((x 2 -ie)m 2 ), (12) 

where m is an arbitrary mass scale which is fixed, for defineteness, on the fermion mass. 
Then, up to an additive constant, 

F=^L— * _ y! + ln [ X '- X Xf_'-% ). (13) 

Arc \ m^\Xf — XfY[Xi — Xi) A [Xf — Xi) z {Xj — Xi) A 1 

Substituting fll3|) into ([□]), we get our final aswer for /^-evolution: 



G f , x j Xi, x^ 



Z(x f - Xi) 2 (x f - Xi 



1 &(0-/*>) 



m 4: (xf — Xf) 2 {xi — Xi) 2 (xj — Xi) 2 {xf — Xi) 2 



x 



G/3 (x f ,Xf,Xi,Xi). (14) 

The normalization Z is infinite before the ultraviolet renormalization. After the renormal- 
ization it is scheme-dependent and calculable order by order in perturbation theory. We will 
not need its value in what follows. 
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3 The Bound State Parameters And The Four-Fermion 
QED Green Function 

The four-fermion QED Green function contains too much information for one who just going 
to calculate bound-sate parameters. Ona can throw away unnessesary information by putting 
senter of mass space-time coordinate of ingoing pair and relative times of both ingoing and 
outgoing pairs to zero: 

G( et )/3(£,x, r', r) = Gp (x/(t,x, r'),x f (t,x, r'), Xj(r'), Xj(r')) , (15) 

where the space-time coordinates depend on a space-time coordinate of the center of mass 
of the outgoing pair (t, x) and a relative space coordinate of outgoing (r') and ingoing (r) 
pair. In the case of equal masses 

z/ = (t,x+~), X/ = (t,X--), 

Xi = (0,^), ^=(0,~). (16) 

G(et)p still contains an unnecessary piece of information — the dependence on the center 
of mass space coordinate. The natural way to remove it is to go over to momentum represen- 
tation and put the center of mass momentum to zero. In coordinate representation, which 
is more convenient for gauge invariance check, we define the propagator Dp of the fermion 
pair: 

G (et) p(t, x, r', r) = Dp(t, r', r)S(x) + . . . , (17) 

where dots denote terms with derivatives of <5(x). It is natural to consider Dp as a time 
dependent kernel of an operator acting on wave-functions of relative coordinate. In what 
follows we will not make difference between a kernel and the corresponding operator. The 
naturalness of the above definition of the propagator is apparent in the nonrelativistic ap- 
proximation: 

e i2mt Dp(t) « ]T#(t)e- i£o 'P(£ ), (18) 

Eo 

where the summation runs over the spectrum of nonrelativistic Coulomb problem and P(Eq) 
are the projectors onto corresponding subspaces of the nonrelativistic state space. One can 
obtain ( 18|) keeping leading term in a-expansion of the lhs if one will keep t oc I /a 2 and 



r', r oc 1/a (see || The subscript on E is to denote that it will get radiative corrections 
(see below). The exponential in the lhs is to make a natural shift in energy zero. In what 
follows we will include the energy shift in the definition of Dp(t). 

The next step in calculation of radiative corrections to the energy levels is a crucial one: 
one should make an assumption about the general form of a deformation of the t-dependence 
of the rhs of (0) caused by relativistic corrections. A naturall guess and the one which leads 
to the generally accepted rules of calculation of the relativistic corrections to the energy 
eigenvalues (see, for example (4j) is to suppose that one can contrive oscillating part of 
the exact propagator Dp from the rhs of fll8|) just shifting energy levels and modifying the 
operator coefficiens P(E ): 

Dp{t)= £ 9(t)e~< Eo+AE °) t Pp(E + A Eo ) + ..., (19) 

e +a Eq 
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where dots denote terms which are slowly-varying in time (the natural time-scale here is 
1/Eq). The additional subscript f3 on Pp is to denote that oscillating part of Dp(t) can 
acquire a gauge parameter dependence from relativistic corrections. 

Let us see how one can use eq.(^) in energy level calculations. It is quite sufficient to 
consider Dp(t) on relatively short times when A^ t -C 1, E^t ~ 1. For such times one can 
approximate Dp expanding the rhs of eq. ([T9|) over A Eo t: 

D p {t) « E^e^'E^W. (20) 

E k 



where 



\k 



4\ E o) = E i - J ^ L Pp(Eo + A Bo ). (21) 



k\ 



An extraction of these objects from the perturbation theory is an interim step in the level shift 
calculations. (Here we should mention that in calculational practice Ap(Eo) are exracted 
in momentum representation — i.e. not as coefficients near the powers of time but as the 
ones near the propagator-like singularities (E — E + ie)~( k+1 \) To come nerer to the level 
shift values, useful objects are 

Af=^Af{E Q )i h k\. (22) 

Eq 

Namely, as notations of fl2~T|) suggest, eigenvalues of should be equal to normalizations 
of bound state wave functions which are driven from unit by relativistic corrections while 
the eigenvalues of Ap should be energy shifts to the k-th power times corresponding nor- 
malizations. Thus, the eigenvalues of 



■4 o 



1 A (k) , A (k) 
A j3 A fi 



4(0) 
A f3 



2 

should be just energy shifts to the k-th power. Thus, we define 



(23) 



Sp = 4 0) (24) 

to be the energy shift operator: its eigenvalues are the energy level shifts caused by relativistic 
corrections. Our aim is now to check /3-independence of Sp eigenvalues. 

Some notes are in order: If the conjecture flT9p is true, Ap should commute with 
and the following relation should hold: 

Sf = [Sp] k (25) 

This relation was suggested as a check of the cojecture (fOJ) in and, to our knowlege, has 
never been checked. Another thing to note is that relativistic corrections affects the form 
of the scalar product of wave functions and, thus, one shoud add a definition of operator 
products to the formal expressions (|23| ) , (|25|) . But the level of accuracy to which we will 
operate permits us not to go into this complication and use the operator products as they 
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are in the nonrelativistic approximation — i.e. as the convolution of the corresponding 
kernels. 

The way to the gauge invariance check of the energy shift calculations is clear now: 
Using the gauge evolution relation ( |14D one should find the /3-dependence of Sp and then of 



its eigenvalues. As Sp is defined in 



through An 



(*). 



s which are, in turn, defined in 



( p0|) through the propagator Dp, the first step is to simplify ([14]) to the reduced case of zero 
relative time and total momentum of the fermion pair: 



Dp(t, r',r) 



7(4* 2 )) 



[1 - ((r' 
Z 



m 2 r' 2 m 2 r 2 



v) 2 /(At 2 )) 



x 



Dp (t,r',v). 



(26) 



The factor in the square brackets of the second line is time- independent and futher factorizible 
on factors depending on either ingoing or outgoing pair parameters. This reduce the influence 
of this factor to a change in the normalization of states. Being interested in gauge invariance 
of energy shifts, we omit this factor in what follows. Let us turn to the analysis of the 
influence of the factor in the first line of (^). 

This factor is close to unit in the atomic scale r', r ~ 1/a, t ~ 1/a 2 . We will use its 
approximate form: 

(27) 



a 



r'r 



Factor ^ 1 + —(/?- /3 ) — + 0(a 



2tt v " t 2 



(k) 

One can read the dependence of A \ on (3 from 



A {k) 



4?- 



a 



(P ~ AO 



2vr(fc + l)(fc + 2 



■tA 



(k+2) 
Po 1 



(28) 



where r is the vector operator of relative position of interacting particles. The mixing of 
different Ap's with a change in the gauge parameter is due to the presence of 1/t 2 in the 



(29) 



rhs of (|27|). Finally, using the definition (gj), relations (|25D and the fact that 

A® w 1 



in the nonrelativistic approximation one can derive the following /5-dependence of Sp\ 



s, 



00 



(30) 



Treating the term in the last line of the rhs of the above relation as a perturbation, one can 
get an approximate value of the /3-dependent piece of the energy shift just averaging the 
perturbation with respect to the corresponding eigenstate of Sp . 

Thus, we get for the leading order of /5-derivative of an energy shift the following repre- 
sentation: 



d_ 



A, 



a 
2^ 



- \ (S L rSlr) - \ (rSlrS L 



(31) 
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where (...) means averaging with respect to the corresponding nonrelativistic eigenstate and 
the subscript L means the leading order in a-expansion. 

Eq.(|3T|) is sufficient to define an order in a in which the energy shifts become gauge 
dependent: 

(! A *)r aU (32) 

Here we have taken into account that r ~ I /a and Sl ~ a 4 . 

To have a gauge dependence in any observable is clearly unacceptable. In the next section 
we will see how one should correct the above procedure of energy shift extraction from the 
QED Green function to get rid of the gauge dependence of energy shifts. 



4 A Way Out 



The procedure recalled in the previous section is based on the conjecture (^9|). A consequence 
of this conjecture is the gauge dependence of energy shifts of fl3~l~|). One can conclude that 
the conjecture is wrong. In particular, as one can infer from eq . (p6|) , the operator coefficients 
near the oscillating exponentials in ([B]) shoud get a time dependence from relativistic cor- 
rections. Even if in some gauge they are time independent, the gauge parameter evolution 
should generate a dependence which in the leading order in a reduce itself to the following 
replacement in (|T9|): 

Pp(E + A Eo ) -> P P {E + A Eo ) + 5^1. (33) 

That S /3 (_E , ) has nothing to do with energy shifts but will give contributions to A^\E o ys 
from eq. (JZ0|) . Being gauge dependent these contributions lead to the gauge dependence of 
energy shifts. 

The way to the correct procedure is to through away terms like H^Ef^/t 2 prior to the 
definition of the energy shift operator. Thus, a necessary step in the process of extracting 
energy shifts from the QED Green function (and the one which necessity is not recognized 
in the stanard procedure) is to calculate and subtract contributions like the last term in the 
rhs of fl33l ) from the propagator of the fermion pair. 

Below we report on a calculation of E^-Eq) from fl3"5|). The most economical way to cal- 



culate it is to note that the energy dependence of the Fourier transform of the corresponding 
contribution to the propagator is 

(E-E )ln(-(E-E + is)) (34) 

and that it comes from diagrams describing radiation and subsequent absorption of a soft 
photon with no change in the level Eq of the radiating and absorbing bound state. Similar 
contributions (with another power of energy before the log) are well known for the propagator 



of a charged fermion |TU 



The first step in our calculation is to present the pair propagator in the following form: 
D p (t) « (^^(t, ^) e_ierA(0) ) A=0 . ( 35 ) 
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where L s is the same as in (|3|) except a restriction on the momentum of photon propagator 
- the range of its variation is restricted to the soft region which border is of order of atomic 
binding energies; the exponentials with gauge potential are originated from the ones in ([?]); 
D inv is a descendant of Gi nv from (0): to go over from Gi nv to D inv one should make all 
pairing of non-soft photons in Gi nv and all the reductions of space-time coordinats which was 
involved in going over from the Gp of (|l|) to the Dp of (|17D ; at last, all gauge potentials in 
(R5|) are taken at zero of space coordinate in accord with the <5(x) of eq. (|TTD . The difference 
between the lhs and the rhs of eq.fl35|) does not conribute to the term under the calculation. 

The leading in the nonrelativistic approximation contribution to D inv is the same as 
for Dp — it is just the propagator of the nonrelativistic Coulomb problem. We explicitly 
calculate the leading contribution to the dependence of D inv (t,A) on the gauge potential 
in its expansion over soft momenta of the external photons. Not surprisingly, the dipole 
interaction of the pair with the external photon field arise in this approximation: 

D mv (t, A) « [i^ -H c + er£(t)j , (36) 

where H c is the hamiltonian of the nonrelativistic Coulomb problem and £ is the strength 
of the electric field: 

£(t) = -A(t)+VA (t). (37) 
Substituting (RBI) in (RH) and keeping terms with only one soft photon propagator we get 



expressions which sum contains the term under calculation 

)a=o ' 



e 2 (L s rA{t)D nr (t)rA(0)) A=QJ (3* 



e 2 (L s J d n dr 2 D nr (t - r 1 )r^(r 1 ) J D nr (r 1 - T 2 )r£(r 2 )D nr (r 2 )^ _ , (39) 
ie 2 (L s J dr(D nr (t- r)r£(T)D nr {r)r A(0) - (40) 



rA(t)D nr (t - T)r£(r)D nr (r)) 

/ A=0 

where D nr (t) is the propagator of the nonrelativistic Coulomb problem from the rhs of 

The next step is to pick out a contribution of a level E in (|38D ,(|39D, (|40|) . That is 
achievable by the replacement 

D nr (t) - e- iEot 9(t)P(E ). (41) 



The last ingredient that one needs to calculate (pq),(p9|),(p)|) is the time dependence of 
the soft photon propagators. It can be deduced from (0) as 



2 . .2\ 



2' 

"2; 

5i 



{L s £i(ti)£j{t 2 )) = 9 ((h - t 2 f > if) (42) 



S 



Here the ^-functions are to account for the softness of the participating photons (t c ~ 1/E ). 
Taking (f42"D into account we get the following contributions from (|38D ,(|39D, (|40|) : 

(38) - ^9(t)e- iE °^ (-1 + ^ rP(£ )r, 
© - ^6(t)e- iE ^P(E )vP(Eo)rP(E ), 

& - ^6(t)e- lEot -t(P(E )rP(E )r-rP(E )rP(E )). (43) 

t Z 71 

The sum of the above terms yields the result of our calculation: 

Zp(E ) = -(lp(E )rP(E )rP(E ) + (-l + ^-)rP(E )r + 

i(P(E )rP(E )r - rP(E )rP(E )f) . (44) 

One can explicitly check that (3- dependence of E j g(£ , ) is the right one — i.e. if one 
subtracts the S-term from the propagator before the definition of the energy shift operator, 
the latter becomes gauge independent. Another observation is that the E-term cannot be 
killed by any choice of the gauge (in contrast to the case of charged fermion propagator 
where an analogous term is equal to zero in the Yennie gauge). 

Summing up, in this paper we derived a relation between QED Green functions of different 
gauges. We used it to check the gauge invariance of the energy shift operator. It turns out 
to be gauge dependent. This fact forced us to recognize that energy shifts are not one, and 
the only one, source for the positive powers of time near the oscillating exponentials in the 
propagator of the pair. We found a particular additional source of the positive powers of 
time which is responsible for the gauge dependence of the naive energy shift operator. We 
conclude by an observation that at the moment we have not a clear definition of the energy 
shift operator — to get it one needs a criterion for picking out contributions to the positive 
powers of time originating from the energy shifts. 

The author is grateful to A. Kataev, E. Kuraev, V. Kuzmin, A. Kuznetsov, S. Larin, 
Kh. Nirov, V. Rubakov, D. Son, P. Tinyakov for helpful discussions. This work was supported 
in part by The Fund for Fundamental Research of Russia under grant 94-02-14428. 
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